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Abstract 

The curvature effect on the electronic states of a deformed cylindrical conducting 
surface of variable diameter is theoretically investigated. The quantum confinement 
of electrons normal to the curved surface results in an effective potential energy 
that affects the electronic structures of the system at low energies. This suggests 
the possibility that ballistic transport of electrons in low-dimensional nanostructures 
can be controlled by inducing a local geometric deformation. 
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1 Introduction 



Recent advances in the manipulation of nanostructures have enabled the fab- 
rication of reduced-dimensional quantum systems with novel geometry [1-11]. 
The understanding of their basic properties and the accurate modeling of their 
electronic structures are of vital importance for the manufacture of nanode- 
vices and their applications. From the theoretical viewpoints, the peculiar 
features exhibited by quantum systems confined to low- dimensional curved 
geometry are of interest. These peculiarities arise from an interplay between 
geometry and quantum physics; in fact, when the electron is strongly confined 
to a smoothly curved surface, it experiences an effective potential energy whose 
magnitude depends on the local curvatures along the surface [12-24]. Due to 
this curvature effect, the electrons can not move around freely on the surface 
even in the absence of impurities or other interacting entities. This implies that 
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the quantum transport of low- dimensional nanostructures can be controlled 
by altering the local geometric curvature. 



In the present work, we theoretically investigate the curvature effect on the 
quantum properties of electrons confined to a cylindrical surface. The local 
deformation of the cylindrical surface significantly affected the spatial profile 
of the electron eigenstates in the lowest-energy region. This indicates the oc- 
currence of a curvature-induced alteration in the ballistic electron transport 
in nanoscale cylindrical surfaces. 



2 The confining potential approach 

In quantum mechanics, the motion of quantum particles constrained to a two- 
dimensional curved surface is described by one of the two formalisms given 
below. One is the intrinsic quantization approach, in which the motion is 
constrained to the surface a priori; namely, a classical Hamiltonian is firstly 
constructed from coordinates and momentum intrinsic to the surface, following 
which the system is quantized canonically. The other is the confining potential 
approach, in which the particle is assumed to be confined by a strong force 
that acts normal to the curved surface. In this approach, the quantization of 
the motion perpendicular to the curved surface results in an effective poten- 
tial that depends on the local surface curvature. Among the two formalisms, 
we employ the latter one since it offers a physically more realistic model of 
quantum confinement to curved surfaces. (In fact, in any real physical system, 
constrained motion is the result of a strong confining force.) It is mentioned 
that the confining potential approach was initially suggested by da Costa [13] 
; it has been successfully applied to quantum mechanical problems involving 
novel geometries [25-30]. 

Let (a; 1 , x 2 , a; 3 ) be a three-dimensional curvilinear coordinate. This allows the 
parameterization of a curved surface of interest by r = r{x 1 ,x 2 ). Then, a 
point p in this space can be determined through the relation 



where n(x 1 ,a; 2 ) is the vector normal to the curved surface. We now introduce 
a confining potential V(x 1 , x 2 , x 3 ). After a proper limiting procedure [13], the 
potential becomes V = if x 3 = and V = oo otherwise. This allows us to 
separate the x 3 dependence in the Hamiltonian of the confined system and 
eventually provides the Schrodinger equation for curved surfaces: 
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Fig. 1. Schematic illustration of a deformed cylindrical surface. The parameter a de- 
termines the spatial extent of the deformed part along the z-axis, and j3 determines 
the extent of a bulge {(3 > 0) or constriction (/? < 0). 

where, 
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In Eq. (2), a(x l , x 2 ) is the wave function of the confined particles, and H{x l , x 2 ) 
and D(x\x 2 ) are the mean and Gaussian curvatures of the surface, respec- 
tively [31]. The occurrence of the non-trivial potential term, —(H 2 — D) in 
Eq. (2), is a direct consequence of the quantization of the motion normal to 
the surface. It is emphasized that this potential term depends only on the 
local surface geometry and not on the mass or charge of the particle. 



3 Model and methods 



3. 1 Schrddinger equations for deformed cylindrical surfaces 



Figure 1 gives a schematic illustration of a cylindrical surface subject to local 
deformation. 1 This curved surface is parameterized by 



r = r(z, 4>) = [R(z) cos 0, R(z) sin 0, z\ , 
where the radius of the cylinder is assumed to vary with z as 

-2z 2 ' 



R(z) = R 



1 + (3 exp 



® 2R o, 



(4) 



(5) 



1 Similar cylindrical surfaces have been considered in Refs. [25], [28] and [32]. 
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Hence, the geometry of the surface in question is determined by the two pa- 
rameters: the parameter a determines the spatial extent of the deformed part 
along the z-axis, and (5 determines the extent of a bulge {(5 > 0) or constriction 

(P < o). 

Our current objective is to deduce the electron eigenstates of the curved sur- 
faces introduced above. It is noteworthy that, due to the rotational symmetry 
of the surfaces, the Schrodinger equation (2) can be further simplified by means 
of the variable-separation method. This is seen by substituting 



a{z,<j>) 



7]{z) 



im(f> 



'R(z)f(z) V27T 
into Eq. (2), with the definition 
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yjl + (dR/dz) 2 

Then, we obtain the reduced Schrodinger equation for rj(z) as 
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(11) 
(12) 



Consequently, the problem is reduced to solving the differential equatioin (8) 
with respect to rj(z). By observing the a- and (3- dependences of r](z) , we can 
clarify the curvature effect on the electronic structures of deformed cylindrical 
surfaces. 
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3.2 Tight-binding approximation 



In the actual calculations in Eq. (8), we have employed the tight-binding 
approximation; i.e., the continuous variable z is discretized into a set of discrete 
numbers i (1 < i < N) with an equiseparation a. The resulting difference 
equation is 

to( 2 fi ~ v i)Vi ~ toffiVi+i + Vi-i) = £Vi, (13) 



where the constant to = h 2 /(2m*a 2 ) denotes the hopping energy between 
neighboring sites in an undeformed region. Nonzero surface curvature at the 
deformed region manifest itself in spatial modulation of fa and Vi given in 
Eq. (13). In fact, /j = and = const, if the cylinder in question is flat (i.e., 
R = const.). 

The differential equation (13) applies to real deformed cylindrical nanosurfaces 
in which the surface density of constituents (atoms or molecules) is homoge- 
neous over the whole surface. In other words, the separation of neighboring 
constituents needs to assume nearly constant so that the hopping energy re- 
mains to be constant. This condition, however, may violate when an external 
deformation against the original flat cylindrical surface is magnified. In the 
latter case, deformation-induced change in can be taken into account by 
rewriting Eq. (13) as 

*o(2/i 2 - Vi)r}i - t ii+1 ffrj i+1 - tu-iffrji-! = erji. (14) 



Here, t^ represents the spatially dependent hopping energy, which decreases 
exponentially with the separation as 



Uj = t e~^' a \ (15) 
where 



nj = J 9u(z)dz, 0n (z) = yj-p. (16) 

Obviously, if = a for all i and j, t^ = to so that Eq. (14) reduces to Eq. 
(13). In both the conditions mentioned above, the radius Rq of an undeformed 
part of the cylindrical surface is set to R = 50 and the length L of the cylinder 
is L = 1000 in units of a. 
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Fig. 2. Contour plots of the localization length £ of the lowest-energy eigenstate 
with (a)-(b) m = and (c)-(d) m = 1. In the two left panels, the hopping energy 
to is assumed to be constant, while in the two right panels, the spatial variation of 
Uj is taken into account (see Eq. (15)). The eigenstate is strongly localized in the 
dark region and extends over the entire system in the light region. 



4 Results and discussions 



4-1 Localization length 



In order to clarify the curvature effect on the spatial profile of \r](z)\ 2 , we have 
calculated the localization length £ of the lowest-energy eigenstate denned by 
[33] 



v 



i=i 



(17) 
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This quantity provides a measure of the spatial extent of the wavefunction 
in question. 2 Hence, by examining the dependence of £ on the geometric 
parameters a and (3, we obtain an understanding of the geometric effect on 
the electron system on cylindrical surfaces. 

Figures 2(a)-(d) show the contour plots of the localization length £ of the 
lowest-energy eigenstate with (a)-(b) m = and (c)-(d) m — 1. In the two left 
panels, the hopping energy t is assumed to be constant (see Eq. (13)), while 
in the two right panels, the spatial variation of Uj is taken into account (see 
Eq. (14)). Henceforth, we refer to the two models described by Eqs. (13) and 
(14) as models A and model B, respectively. 

We see from the four panels that the ground-state eigenfunction may be 
strongly bounded (dark region) or may extend over the entire system (light 
region) depending on the values of a and (3. We first consider the result of 
model A; it is evident that in both the cases of m = (Fig. 2(a)) and m — 1 
(Fig. 2(c)), the a dependence of £ is completely different in the regions (3 > 
(bulging deformation) and (3 < (constricted deformation). With m — 0, for 
instance, £ for (3 > tends to increase with a and eventually becomes equal to 
the system length L at a ~ 1.5. This indicates that the ground-state becomes 
extended (£ ~ L) when reforming the bulging part to be more stretched in the 
direction of the z-axis. On the other hand, £ for f3 < is almost invariant with 
regard to changes in a, and its value is always smaller than L. As summarized, 

(i) the crossover behaviour of £ from the bounded to extended states is ob- 
served only in the case of bulging deformation ((3 > 0), and within this range, 

(ii) the magnitude of £ for a fixed (3 rapidly increases at a ~ 1.5. These two 
findings suggest the possibility of controlling the ballistic electron transport 
in cylindrical nanostructures by introducing a subtle geometric deformation. 

A similar crossover behaviour of £ was observed in the case of m = 1 (Fig. 2(c)). 
Nevertheless, in contrast to the case of m — 0, the crossover occurs for neg- 
ative (3, not for positive. In fact, we have confirmed that when m > 1, the 
extended low-energy states occur only at f3 < 0, which is in contrast to the 
case of m = 0, where the extended states occur only at (3 > 0. This difference 
as well as the global behavior of £(a, (3) depicted in Figs. 2(a) and 2(c) can be 
accounted for by considering the contribution from each term in the expression 
for v(z) (see Eq. (9)); detailed analyses will be presented elsewhere [35,36]. 

While the localization length £ in the model A is responsive to the variations 
of a and (3, that in model B exhibits only a slight dependence on the values of 
a and (3. This is illustrated in Fig. 2(b) for m = and Fig. 2(d) for m = 1. In 
both cases, £ decreases rather slowly with an increase a and (3, which implies 
that a geometric deformation has almost no contribution to determining the 

2 Instead of Eq. (17), there is an alternative definition of localization length asso- 
ciated with the Lyapnov exponent; see Ref [34] for example. 
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Fig. 3. The upper panels: Spatial profiles of |j?(^)| 2 of the ground-state with m = 
for model A (the solid line) and for model B (the dotted line). The middle (lower) 
panels: The z-dependence of the off-diagonal (diagonal) element of the Hamilto- 
nian matrix for model A (solid) and B (dotted). The parameters are set to be 
(a,/3) = (0.8, 1.0) in the left panels, and (a,/3) = (1.7, 1.0) in the right panels. 

spatial profiles of the eigenstates. 



4-2 Spatial profile of the wavefunction 



To understand the difference in the relevance of geometric deformation to the 
behavior of £ between in models A and B, we extract the spatial profile of 
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the square amplitude of the ground-state eigenfunction for m = by fixing 
a and (3 to be certain values. Figures 3(a) and 3(b) are plots of the profile 
of |^(^)| 2 for model A (the solid line) and for model B (the dotted line) with 
different values of (a, (3). The parameters are set as (a, (3) = (0.8,1.0) in 
the left panels and (a, (3) = (1.7, 1.0) in the right panels. We have also plotted 
the ^-dependence of the diagonal and off-diagonal elements of the Hamiltonian 
matrix for the two models in Figs. 3(c)-(f), where to is taken as units of energy. 

We see from Figs. 3(a) and 3(b) that the wavefunction for model A is extended 
at a = 0.8, whereas it is spatially bounded around the deformed region (z ~ 0) 
at a = 1.7. Indeed, this behavior is completely consistent with the contour 
plot of £ depicted in Fig. 2(a). In contrast, the wavefunction corresponding 
to model B exhibits almost no variance in its spatial profile for a different 
value of a; it is fairly extended along the z-axis with a low (almost negligible) 
amplitude around the deformed region. 

The disappearance of the bound state in model B is qualitatively understood 
by observing the z-dependence of the off-diagonal elements tijff, as shown 
in Figs. 3(c) and 3(d). In both plots, double- well structures whose bottoms 
are located at z — rkz w appear. It is noteworthy that the magnitude of the 
double well of model B is larger than that of model A. This is because in the 
model B, the value of the hopping energy tij at z — ±z w becomes smaller than 
to, since the separation becomes larger than a (see Eq. (15)). Then, the 
quantum hopping across point ±z w is slightly weak; thus, the wavefunction in 
model B tends to have a finite amplitude only within the undeformed region 
(\z\ > z w ) in order to make the energy lower. As a consequence, the geometric 
deformation in model B neither engenders spatially bounded eigenstates as 
done in model A nor contributes significantly to the spatial profile of the 
lowest-energy eigenstates. This scenario indicates that the low-energy bound 
states observed in model A are difficult to realize by inducing an external 
mechanical deformation to a flat cylindrical conducting surface. We emphasize 
that the abovementioned findings regarding model B serve as a complementary 
result to the existing studies of the geometric effect on the electron transport 
in cylindrical surfaces, where the geometric conditions to realize the bound 
states were considered [25,28]. 



5 Conclusions 

We have investigated the effect of local geometric deformation on the electronic 
states that are strongly confined to a thin cylindrical surface. The spatial pro- 
file of the square amplitude of the lowest-energy eigenstates were numerically 
evaluated by the tight-binding approach; this was followed by the extraction 
the localization length of the eigenstate as a function of the geometric param- 
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eters a and (3. The results indicate that, under certain conditions, a subtle 
geometric deformation could induce a drastic change in the ballistic electron 
transport along nanoscale cylindrical surfaces. We hope that our findings prove 
to be a fundamental basis for the development of quantum devices based on 
low-dimensional nanostructures. 
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